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ABSTRACT 

We present 13 cm and 20 cm radio observations of the magnetic chemically peculiar star 
CU Virginis taken with the Australia Telescope Compact Array. We detect two circu- 
larly polarised radio pulses every rotation period which confirm previous detections. In 
the first pulse, the lower frequency emission arrives before the higher frequency emis- 
sion and the ordering reverses in the second pulse. In order to explain the frequency 
dependence of the time between the two pulses, we construct a geometric model of the 
magnetosphere of CU Virginis, and consider various emission angles relative to the 
magnetic field lines. A simple electron cyclotron maser emission model, in which the 
emission is perpendicular to the magnetic field lines, is not consistent with our data. 
A model in which the emission is refracted through cold plasma in the magnetosphere 
is shown to have the correct pulse arrival time frequency dependence. 

Key words: radio continuum: stars — stars: individual: CU Virginis, magnetic fields, 
rotation — radiation mechanisms: non-thermal 



1 INTRODUCTION 

CU Virginis (HD124224, hereafter CU Vir) is one of the 
best studied Ap stars. It is very nearby (distance of only 
80 pc) and is a fast rotator (period of around 0.52 days) 
(|Deutschl [l952h . Its strong magnetic field (polar magnetic 
field of ^3 kG) places it in the category of Magnetic Chem- 
ically Peculiar (MCP) stars. Like other MCP stars, its 
effective magnetic field has been observed to vary with 
rotational phase (Bor ra &; Landstreetl Il980f). In addition , 
MCP stars are known radio emitte rs (|Linskv et al.l tl992: 
iLeone et all 1 19941 : [Drake et al.lll987h and their radio emis- 
sion has been observed to modulate with rotational phase 
(Leone & Uma na 1993: iLeond Il99lh . This quiescent radio 
emission is believed to be gyrosynchrotron radiation emit- 
ted by mildly relativistic (Lorentz fac tor of 7 < 2) elec - 
trons trapped in t he magnetosphere ([Drake et al.l Il987h . 
ITrigilio et a 1. (2004) constructed a three-dimensional model 
to explain the rotational modul ation of the quiesc ent radio 
emission from a MCP star and iLeto et aD (2006) success- 
fully applied it to CU Vir. 

Circular polari sation in the radio emission was detected 
at the 10% level bv lLeone et al.1 |l996) which is expected for 
the gyrosynchrotron emission. Subsequently, ITrigilio et al.l 



(2000) discovered that CU Vir produces two 100% circularly 
polarised radio pulses every rotation period which cannot 
arise from the gyrosynchrotron emission mechanism. 

The observed radio flux for CU Vir across one rota- 
tion period is hereafter described as the "pulse profile" . The 
two components in the pulse profile are referred to as the 
"leading" and "trailing" pulses. Similar pulse profiles have 
been obser ved from CU Vir in the 13 cm and 20 cm observ - 
ing bands (jTrigilio et al.ll2000U2008ll201ll : iRavi et al.ll2010h 
where observations were separated by almost a decade. CU 
Vir is the only MCP star that has been observed to emit 
radio pulses. 

ITrigilio" et al.l (|2000h proposed electron cyclotron maser 
(ECM) as the emission mechanism for the pulsed radio emis- 
sion. Pl^sjna^a^ation due to Langmuir waves was ruled 
out (iTrigilio et al.l l2008h because of the low electron den - 
sity in the emission region as modelled bv lLeto et all (2006) 
and the observation that requires the radiation to be tightly 
beamed. ECM can produce radiation with high directivity 
and 100 % circular polarizatio n. In the "loss cone" model of 
ECM (jMelrose fe Durklll982h . ECM emission forms a hol- 
low cone at a large angle (70° to 85°) to the magnetic field 
lines, and occupies a narrow frequency band close to a har- 
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monic of the cyclotron frequency. Although ECM is narrow 
band, the relatively wide band of the observation can be 
explained if the emission region spans a range of magnetic 
field strengths. 

A simple ECM model does not explain all the features of 
the pulsed radio emission from CU Vir. For instance, in the 
20 cm observing band, the leading and trailing pulses are 
separated by ~ 5 hours. In the 13 cm observ ing band, the 
leadin g pulse was not seen in observa tions by Trigilio et al.l 
(2008), but was seen a decade later bv lRavi et al] (|2010l ). In 
contrast to the variability of the leading pulse, the trailing 
pulse shows remarkable stability over more than a decade 
in both observing bands. The trailing pulse at 13 cm is al- 
ways observed to arr i ve ear lier than the corresponding 20 cm 
pulse. iTrigilio et all ([2008T ) postulated that this has a geo- 
metric explanation, which is related to where the pulse orig- 
i nates in the magnet osphere. The discovery from I Ravi et al.l 
that the leading 13 cm pulse arrives later than the 
corresponding 20 cm pulse, agrees with this interpretation 
and suggests that both pulses come from the same magnetic 
pole of the star. However, until now, existing data do not 
provide enough frequency coverage to study the pulse pro- 
file in detail as a function of observing frequency. Here we 
present new data with more complete frequency coverage 
over the 13 cm and 20 cm observing bands and demonstrate 
with our modelling that an ECM emission model, without 
propagation effect s, cannot exp lain the observations. In a 
recent paper, ITrigilio et alj (|2Q1 lh suggest the frequency de- 
pendence of the pulse arrival time arises from frequency de- 
pendent refraction of the radiation through the stellar mag- 
netosphere. In this paper, we have analysed the frequency 
dependent refraction model in more detail. 

The structure of the paper is as follows. In §2 we present 
the observing and analysis techniques used to detect time 
varying emission from CU Vir with the Australia Telescope 
Compact Array (ATCA) and present our observational re- 
sults. In §3 we discuss our modelling and show how it relates 
to our observations. 



2 OBSERVATIONS AND DATA ANALYSIS 

We observed CU Vir with the ATCA on six separate oc- 
casions (2009 December 23 to 26, 2010 May 19, 2010 June 
15). The observations in December 2009 were only two hours 
each in length because they were taken during unallocated 
telescope time. For these short observations, we predicted 
the arrival time of the trailing pulse using the ephemeris 
from I Ravi et "al . (2010) and timed our observation accord- 
ingly. We were allocated time for the May and June 2010 
epochs and hence were able to observe for about nine hours 
continuously. A summary of the observation parameters is 
given in Table [T] 

All our observations were taken with the Compact Ar- 
ray Broadband Backend (CABB; Wilson 2011) which offers 
a bandwidth of up to 2 GHz per intermediate frequency and 
polarisation. However, the maximum bandwidth in our 20 
and 13 cm observing bands was limited to around 700 MHz 
because our observations were taken before the upgrade was 
fully completed. Furthermore, it was not possible to observe 
at 20 and 13 cm simultaneously. To maximise the apparent 
bandwidth, we switched between the 20 and 13 cm bands 



every few minutes in the 2009 Dec 25 to 26 and 2010 May 
19 epochs. Complex visibility data in four polarisations were 
recorded for 15 baselines with a sampling time of 10 s. 

We also made use of archival (pre-CABB) ATCA data 
to supplement our analysis. Altogether, we analysed eight 
epochs of data at 20 cm and seven epochs at 13 cm as sum- 
marised in Table [1] In column order we provide the date of 
observation, the observing frequencies useqj, the available 
bandwidth, the configuration of the ATCAq the total inte- 
gration time, the phase calibrator used and a reference for 
the data. 

We calibrate d our data using the standard miriad 
software package (Sau lt et al.l [1995)0. The flux amplitude 
scale and the bandpass response were determined from the 
backup ATCA primary calibrator 0823—500. Observations 
of a bright compact radio source (as listed in Table [1]) for 1.5 
minutes every 20 minutes were used to calibrate the com- 
plex gains and leakage between the orthogonal linear feeds 
in each antenna. Radio frequency interference was removed 
from the visibilities using the flagging tool mirflag. 

To obtain light curves, we first imaged the field in the 
standard way. Then, we took the clean components of each 
source, masking the location of CU Vir, and subtracted them 
from the visibility data using the miriad task uvmodel. We 
then shifted the phase centre of the visibilities to the location 
of CU Vir and measured the flux density directly from the 
visibilities by vector averaging their real components in time 
bins of three minutes. 



2.1 Pulse profile 

The 20 and 13 cm Stokes V pulse profiles from our ob- 
servations, along with other archival epochs, are shown in 
Figures [T] and [2] respectively. The Stokes V flux densities 
are plotted agains t rotational phase using the ephemeris 
from ITrigilio et al.l ([201 lh and adjusted for the rotation of 
the Earth around the Sun. We hav e used the ro t ation pe- 
riod, P = 0.52071601 days from ITrigilio et all (|201lh to 
align the light curves and set phase zero, in HJD, to be 
2450966.3601 + EP, where E is the number of epochs. As 
can be seen in Figures [1] and [2] the rotation period is a good 
fit to our data as the pulses are aligned in phase. 

In the May 2010 epoch, we detected both the leading 
and trailing pulses. The time between them is 5.20=b0.05 
hours in the 20 cm observing band. This is similar to the 
pulse separatio n tim e determined by Trigilio et al.l (2000) 
and Ravi et al.l (|2010l ) in earlier epochs. In Figure [3] we have 
juxtaposed measur ement of the effective magnetic field from 
IPyper et al.1 (|l998h with the average Stokes V pulse profiles. 
Assuming a dipolar field structure, and take negative effec- 
tive magnetic field to mean the field lines point into the star, 
the magnetic South pole is closest to our line of sight during 
the shorter of the two pulse separations. 



1 We refer to the 1384 MHz pre-CABB and 1503 MHz CABB 
frequencies as the 20 cm band; and the pre-CABB 2496 MHz and 
CABB 2335 MHz frequencies as the 13 cm band 

2 In all cases a maximum baseline of 
used. The various labels (6A, 6C and 
to different antenna spacings, which are 



6 km was 
6D) refer 
defined at 



http : //www . narrabri . atnf . csiro . au/ operations/ array_conf igurations/ conf ign 
6 http : //www . atnf . csiro . au/computing/sof tware/miriad/ 
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Table 1. Summary of ATCA observation parameters 





Start time 


rrequency 


-Dandwicitn 


Al OA configuration 


• 

Integration time 


Phase calibrator 


Reference 




T TT 
U 1 


^IVlrlZ ) 


(^IVlrlz ) 




(hours) 






1999 May 29 


06:57 


1384, 2496 


128, 128 


6A 


9.3 


B1406-076 


Tridlio et al. (2008) 


1999 Aug 29 


01:12 


1384, 2496 


128, 128 


6D 


9.0 


B1406-076 


Trieilio et al. (2008) 


2008 Oct 30 


21:06 


1384, 2368 


128, 128 


6A 


9.0 


B1416+067 


Ravi et al. (2010) 


2009 Dec 23 


18:13 


1503 


742 


6A 


1.2 


B1406-076 


This work 


2009 Dec 24 


18:11 


2335 


678 


6A 


1.6 


B1406-076 


This work 


2009 Dec 25 


20:01 


1503, 2335 


742, 678 


6A 


1.9 


B1406-076 


This work 


2009 Dec 26 


20:13 


1503, 2335 


742, 678 


6A 


2.4 


B1406-076 


This work 


2010 May 19 


08:02 


1503, 2335 


742, 678 


6C 


8.2 


B1416+067 


This work 


2010 Jun 15 


06:41 


1503 


742 


6C 


9.2 


B1402-012 


This work 



Table 2. Phases of the leading and trailing pulse in the 19 May 
2010 epoch 



Frequency 


Leading pulse 


Trailing pulse 


Pulse separation 


(MHz) 


arrival phase 


arrival phase 


(hours) 


1.375 


0.360 ± 0.002 


0.782 ± 0.002 


5.27 ±0.05 


1.675 


0.366 ± 0.002 


0.778 ±0.002 


5.15 ±0.05 


1.875 


0.372 ± 0.006 


0.776 ±0.006 


5.1 ±0.1 


2.275 




0.76 ±0.02 




2.475 




0.76 ±0.02 




2.675 




0.756 ±0.006 





We did not detect the 13 cm leading pulse, which had 
been detected in 2008. Although the trailing pulse was de- 
tected in all epochs for which we have data at the relevant 
phase, the peak flux density varied between 6-12 mJy in the 
20 cm observing band, and 7-13 mJy in the 13 cm observ- 
ing band. With our data we can constrain the frequency at 
which the pulse disappears to be between 1.9 and 2.3 GHz 
in the 19 May 2010 epoch. 

Figure [4] shows the pulse profiles, in 200 MHz channels, 
for the 19 May 2010 epoch. The pulse arrival phase is de- 
pendent upon the observing frequency. This dependence is 
clearly observed for the trailing pulse in the 26 Dec 2009 
epoch (Figure [5}. This frequency dependence is stable over 
at least six months since we observe similar phenomenon in 
the trailing pulse in the 26 Dec 2009 epoch. We determined 
the arrival phase of the pulses by cross correlating the single 
epoch pulse and the average pulse, and finding the phase 
of the peak of the cross correlation. The error of the pulse 
arrival phase is the variation of the parameter required to 
increase the reduced x 2 between the single epoch pulse and 
the average pulse by one unit. Table[2]summarises the arrival 
phase for the 19 May 2010 epoch. 




Average 



30th Oct 2008 




23rd Dec 2009 



25th Dec 2009 



26th Dec 2009 



19th May 2010 



3 MODELLING 

Models for the origin o f the pulsed emissio n from CU Vir 
have b een proposed by Trigilio et al.l (|2000l ) , ITrigilio et al.l 
(2008) an d ITrigilio et all (|201 ll ). Weperform numerical sim- 
ulation of each of these models to produce a pulse profile, 
and compare these with our observations. 



Figure 1. 20 cm Stokes V pulse profiles from all epochs listed 
in Ta ble [T] aligned using the rotation period from ITrigilio et al.l 
(l201lh. 
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Figure 2. 13 cm Stokes V pulse profiles from all epochs listed 
m Table □ aligned using the rotation period from ITrigilio et al.l 
(120111 ) 



3.1 The magnetosphere of CU Vir 

The emission models described here are based on astel- 
lar magnetosphere model proposed by [Andre et al.1 (|l988h . 
shown in Figure [6] The cold torus, containin g trapped ma- 
terial from the stellar wind, is an extension by ITrigilio et al.l 

(HqoI). 

In this model, electrons are accelerated by magnetic re- 
connection in the magnetic equatorial plane of the star, and 
propagate along field lines towards the magnetic poles. They 
are reflected through magnetic mirroring; however, electrons 
with sufficiently small pitch angles will instead intersect the 
surface of the star. This results in a loss-cone anisotropy 
in the pitch angles of this population of electrons, allow- 
ing them to produce radio emission through t he electron 
cyclotron maser mechanism (|Trigilio et al.ll2000[ ). 

ECM emission occurs close to a harmonic of the cy- 
clotron frequency: v — sis cy c, where s is the harmonic num- 
ber. We assume the emission to be in the second harmonic, 
s — 2, because the fundamental, s = 1, is generally sup- 
pressed by gyromagnetic absorption of the thermal plasma 




0.2 0.4 0.6 0.8 1 

Phase 



Figure 3. Top: Effective magnetic field B e ff from IPyper et al.l 
(1998) vs. rotation phase. A negative B e ^ means the field lines 
are directed into the star. Middle and bottom: Average Stokes 
V pulse profiles at 13 cm and 20 cm. If we assume CU Vir has a 
dipolar field structure, then the magnetic South pole is closest to 
our line of sight during the shorter of the two pulse separations. 

in weaker magnetic fields as the radiation escapes, and the 
growth rates for the highe r harmonics are too l ow for the 
intensity to be significant ([Melrose &; Dulklll982h . 

The sense of the polarised emission (positive Stokes V; 
right circular polarisation) indicates that it originates en- 
tirely from the north magnetic pole of the star. This may 
be explained by the presence of a quadrupole comp onent in 
the magnetic field of CU Vir ([Trigilio et al.ll200oh . Conse- 
quently, we consider only the north magnetic hemisphere in 
our modelling. 

3.2 Model parameters 

Table [3] lists the parameters we use for our modelling. 
The obliquit y of t he magnetic axis was calculated by 
ITrigilio et aD ([2000h from measurements of the variation 
in effective magnetic field strength given an inclination i. 
iLeto et al . (2006) derived the Alfven radius and thickness 
of the magnetosphere by fitting multifrequency radio flux 
densities to a three-dimensional magnetospheric model sim- 
ulating the quiescent gyrosynchrotron emission. The dipole 
moment is calculated from a value of 3kG for the mag - 
netic fi e ld str ength at both poles from ITrigilio et al.l (|2000l ) . 
lHatzesI (|l997h suggest an off-centre dipole model, with differ- 
ent surface magnetic field strengths at each pole, but these 
values give a similar dipole moment0. We use only the cen- 
tral value for each parameter, ignoring the uncertainty. 

3.3 Emission from field lines parallel to the 
magnetic axis 

The first model for the pulsed emission from CU Vir is from 
Trigilio et al.1 (|2000l l. The emission is assumed to be directed 
as a hollow cone, making an angle of 85° with the local mag- 
netic field. At this point, it was believed that the population 

4 Since the star is treated here simply as a rotating dipolar field, 
the actual position of its surface is not required for the modelling. 
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Figure 6. Cross-section of the magnetosphere of CU Vir used in our modelling (not to scale). The magnetic field is assumed to be 
dipolar within the Alfven surface (heavy dashes), with equatorial radius r a- Beyond this radius, field lines (light dashes) are distorted by 
the stellar wind. Electrons are accelerated in equatorial current sheets beyond the Alfven radius, and travel along field lines within the 
middle magnetosphere (shaded), which divides the inner and outer magnet ospheres. The parameter I describes the extent of the middle 
magnetosphere, indicating the position of the projection (dotted) of the field line marking its outer boundary within the Alfven surface. 
Cooled gas from the stellar wind accumulates in a torus (shaded) near the star. 



Table 3. Adopted parameters for CU Vir. 



Parameter 


Symbol 


Value 


Rotation axis inclination [de- 


i 


43±7 a 


grees] 






Magnetic axis obliquity [de- 


(3 


74±3 a 


grees] 






Stellar radius [Rq] 


R* 


2.2 b 


Equatorial Alfven radius [R*] 


r A 


15±3 C 


Thickness of middle magneto- 


I 


1.5 C 


sphere [R*] 






Dipole moment [Am 2 ] 


m 


5.4 x 10 33 d 



a iTrigilio etall (120001 ) 
b lNorthl (119981) 



c iLeto et al.l (|20Qi 

d see text 

of high-energy electrons would approach the star along field 
lines almost parallel to the magnetic axis, as shown in Figure 
H (left panel). 

Our simulation procedure is as follows: 

(i) Select a rotation phase cj). 

(ii) Determine the magnetic latitude Oe oi emission in the 
direction of the Earth, using: 



sin Oe — sin /3 sin i cos cj) + cos f3 cos i 



(i) 



See appendix lAll for details. 

(iii) Determine the strength of emission in this direction. 
The emission pattern is a hollow cone making an angle of 
85° with the north magnetic axis; i.e., at a magnetic latitude 
of +5°. We assume the emission to have a narrow Gaussian 
profile centred on this latitude. 

(iv) Repeat steps (ij]|(iii) for different rotation phases to 
form a pulse profile. 

A phase of cj) = under the definition used in equation 
[1] corresponds to the point in the rotation of CU Vir when 
its north magnetic pole is most closely directed towards us. 
Therefore, we shift the phase by 0.08 to match the ephemeris 
as described in section [2] We do not attempt to predict the 
absolute flux density, but instead scale the peak of the sim- 
ulated pulse profile to match our observations. 

The result of this simulation is shown in Figure [3 (right 
panel). This model correctly predicts the arrival times of the 
pulses. However, it does not fit the width of the pulses, nor 
their frequency dependence (which had not been observed 
when this model was developed). 

3.4 Emission from dipolar magnetic field lines 

This model for the pulsed emission from CU Vir is from 
Trigi lio et al.l (|2008h . It incorporates results from ILeto et all 
( 2006), who determined parameters of the magnetosphere, 
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Rotation phase 



Figure 7. Schematic diagram (left), and pulse profile of emission (right) for the model of ECM emission from field lines (dashed) parallel 
to the magnetic axis (section 13. 3|) . The simulated emission in the 20 cm (dark arrows) and 13 cm (light arrows) bands arises from different 
locations, but is emitted in the same directions, hence producing identical pulse arrival times (top right), in contrast to the averaged 
observed pulse profile (bottom right). Dotted and solid lines represent the 13cm and 20cm bands respectively. 






2.675GHz — i — . 








2.475GHz — i — . 








2.275GHz — i — 








1.875GHz — i — 




1.675GHz — i — 












1.375GHz — — 











Figure 4. Stokes V pulse profile from the 19 May 2010 epoch. 
Each panel represents the average of 200 MHz bandwidth. The 
noise is higher for the three pulse profiles in the 13 cm observing 
band compared to 20 cm observing band because more channels 
were removed due to RFI. 



including the Alfven radius r a • These values imply that the 
field lines of the middle magnetosphere near the star are 
not parallel to the magnetic axis. We assume them to have 
a dipolar configuration, as shown in Figures \6\ and [9] (left 
panel) . 

Our simulation procedure is as follows: 



Figure 5. Stokes V pulse profile from the 26 Dec 2009 epoch. It 
is clear that the 1.4 GHz pulse lags the 2.7 GHz pulse by ~ 0.03 
phase. 



(i) Select a rotation phase and convert to magnetic lati- 
tude, as in section [331 

(ii) Select a point on the magnetic equator at radius ta- 

(iii) From this point, trace along the field line towards the 
north magnetic pole, following the path that would be taken 
by an electron in the middle magnetosphere. 

(iv) At each point along the field line, determine the cy- 
clotron frequency. When the cyclotron frequency is equal to 
half of the required emission frequency (corresponding to 
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Figure 8. Emission in a hollow cone (dotted) centred on a mag- 
netic field vector ('B-field') at a magnetic latitude of Ob- The 
half-opening-angle of the cone is 5. A sample emission vector is 
shown as a dashed arrow. This vector may be parameterised with 
the angle uj, which describes its position on the cone. Equation 
[2] describes the relationship between uj and the magnetic latitude 
of the emission, Oe- 



emission at the second harmonic), stop and note the mag- 
netic latitude of the direction of the magnetic field at 
this point. 

(v) The emission is assumed to be directed as a hollow 
cone centred on an axis with magnetic latitude Ob- We de- 
termine the parameter uj which describes the position of the 
emission vector on this cone, as shown in Figure [8] 



sin Oe — cos S sin b 



(2) 



sin S cos Ob 

(vi) If there was no solution for uj, the emission power P 
in this direction is zero. Otherwise, the emission power per 
solid angle Q is: 

dP 1 



sin uj cos Ob sin 

This formula is obtained by assuming the emission power to 
be equally distributed around the cone. See appendix I A2I for 
details. 

(vii) Repeat steps [(ij|(vi)| for different rotation phases to 
form a pulse profile. 

(viii) Repeat steps (I}|(vii) for different values of 5. Sum 
the resulting profiles, weighting them with a narrow Gaus- 
sian function centred on 5 — 85°. 

(ix) Repeat steps (IJ|(viii)' for emission frequencies across 
the observation band. Sum the resulting profiles, weighting 
them equally. 

(x) Repeat steps (ij|(ix) for starting points on the mag- 
netic equator with radii between ta and ta~\-1- Sum the 
resulting profiles, weighting them equally. 

The phase and height of the pulse profile are then ad- 
justed as in section [331 The Gaussian function used in step 
(viii) is very narrow (width < 0.2°), intended only to smooth 



over the step size in other parameters, so the width of the 
resulting profile is due to other aspects of the model. 

Results are shown in Figure [9] (right panel) . The simu- 
lated pulses are too wide to fit the observations, and display 



a distinct double-peaked structure which is absent in the real 
data. There is some frequency dependence, but it affects the 
pulse width (narrower at higher frequency) rather than the 
arrival time. 



3.5 Refracted equatorial emission 

This model for the pulsed emission from CU Vir is from 
ITrigilio et a 

D 4201 ih . Rather than a hollow cone centred on 
the local field line, the ECM emission is assumed to be di- 
rected perpendicular to the local field line, and only in the 
two directions which are also parallel to the magnetic equa- 
torial plane. For each longitude in the magnetic equatorial 
plane, this results in emission being visible from two points 
in the emission region. 

The key point of this model is that it includes the 
effect of refraction as the radiation crosses the boundary 
of the cold torus region (Figure [6]). We neglect the effects 
of refraction from density variations within the cold torus, 
or from the interface where the radiation exits the torus. 
The refractive index outside the cold torus is assumed to 
be unity. Inside the torus, the refractive index n c t depends 
on the orientation of the magnetic field relative to the ra- 
diation (|Gurnett &; Bhattacharieell2005h , but if we assume 
them to be parallel it is: 



n ct 



(4) 



For a plasm a density i n the cold torus of 10 9 cm -3 , 
compatible with iLeto et al.l (|2006r ) , the plasma frequency is 
Up = 280 MHz. If we assume that the magnetic field strength 
at the point of refraction is approximately the same as at the 
point of emission, we have the cyclotron frequency v C yc — 
z//2 (if the emission is at the second harmonic). 

If the angle of incidence on with the boundary of 
the torus (see Fi gure I1Q|) is assumed to be 60°, as by 
ITrigilio et al.l (|201lh , the angle of refraction qsr can be found 
with Snell's law: 



smojK — smo^ 

n ct 



(5) 



For n ct < sinoj^, this implies that the radiation is en- 
tirely reflected at the interface, and none is transmitted. 
Under the conditions assumed here, this occurs for v < 790 
MHz, so the pulse spectrum should be truncated below this 
frequency. However, this is below the frequency range of our 
observations, so we cannot test this here. 

Our simulation procedure is as follows: 

(i) Select a rotation phase and convert to magnetic lati- 
tude, as in section [331 

(ii) Convert from the magnetic latitude of the emission 
to the angle of refraction. From our assumption about the 
angle of incidence, this is done by adding 60°. 

(iii) Calculate the refractive index at this frequency with 
equation [4] 

(iv) Determine the angle of incidence with equation [5] 

(v) Convert from the angle of incidence to the original 
magnetic latitude of the emission. From our assumption 
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Figure 9. Schematic diagram (left), and pulse profile of emission (right) for the model of ECM emission from field lines diverging from 
the magnetic axis (section 13.41 The regions in the magnetosphere where emission occurs are marked with solid rings, and the dark and 
light arrows correspond to emission in the 20 cm and 13 cm bands. The simulated emission (top right) is spread over a very broad range 
of magnetic latitudes, resulting in an extended pulse profile quite distinct from our averaged observations (bottom right). Dotted and 
solid lines represent the 13 cm and 20 cm bands respectively. 



about the angle of incidence, this is done by subtracting 
60°. 

(vi) Determine the strength of emission in this direction. 
We assume the emission to have a narrow Gaussian profile 
centred on the magnetic equator. 

(vii) Repeat steps (i)]|(vi) for different rotation phases to 
form a pulse profile. 

(viii) Repeat steps [(i}| (vii) | for emission frequencies across 
the observation band. Sum the resulting profiles, weighting 
them equally. 

The phase and height of the pulse profile are then ad- 
justed as in section [331 The Gaussian function used in step 
is very narrow (width < 0.2°), so the width of the re- 



sulting profile is due to variation with frequency across the 
observation band. 

Results are shown in Figure [10] (right panel) . The pulse 
arrival times and their frequency dependence both fit the 
observations. The observed pulses are wider than the simu- 
lated pulses; this may be due to inherent width of the emis- 
sion beam, variation of the angle of incidence, or further 
turbulent refraction in the cold torus. 



4 CONCLUSIONS 

We observed CU Vir with the ATCA over six epochs and, 
in the 20 cm observing band, detected 100% circularly po- 
larised pulses twice in a rotation period. This confirms that 
the 20 cm pulses are stable over a period of more than a 
decade. We did not detect the 13 cm l e ading pulse that was 
observed a year earlier by iRavi et al. (2010). For our data, 
the leading pulse has a spectrum which cuts off between 
1.9 and 2.3 GHz. We can see a clear frequency dependency 
in the time between leading and trailing pulses across our 
entire observing band. 

We h ave simulated i n deta il the emission models pro- 
posed bv ITrigilio et all (|20QQh . iTrigilio et all (|2QQ8h and 
ITrigilio et al.l (|2Q1 lh . We show that the first two models do 
not fit our observations, the first being excluded by its lack of 



frequency dependence and the second by its grossly dissimi- 
lar pulse profile. We confirm that the third model, involving 
refraction of the emission as it propagates through material 
in the stellar magnetosphere, yields the correct frequency 
dependence and pulse arrival times to fit our observations. 
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U = E-(E-S)S 
V = B - (B • S)S 



with lengths: 



\U\ = sini 
\V\ = sin/3 

The angle between these vectors is 0, which we can then 
find as: 



U V 



\u\\v\ 

E • B - (E • S)(S • B) 
sin i sin f3 



(Al) 



As E, S and B are unit vectors, we can evaluate these 
dot products as trigonometric identities: 



smOs — cos i cos /3 



sin i sin /3 



(A2) 



Rearranging, we obtain: 



sin Be — sin f3 sin i cos + cos /3 cos i (A3) 
which is the same as equation [1] 



APPENDIX A: ANGULAR RELATIONS 
Al Rotation of CU Vir 

Equation [T] connects the rotation phase of CU Vir to the 
magnetic latitude of emission in the direction of the Earth. 
The relation between these is shown, with associated angles 
and unit vectors, in Figure IA1I To derive the relationship, 
we construct vectors: 



A2 Conical emission 



Step (v) of the procedure in section \3A\ assumes a relation, 



equation [21 between the magnetic latitude of emission Be 
and the position uj of that emission on a cone, as shown in 
Figure [8] To demonstrate this, we note that Figure IA1I is 
equivalent to Figure [8] with the following substitutions: 
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Magnetic axis 
M 




Figure Al. Rotation of CU Vir. The magnetic axis M and the 
rotation axis R are inclined by (3. In the reference frame of the 
star, the vector towards the Earth E describes a circle around 
point C on the rotation axis, separated from the rotation axis 
by i, and parameterised by the rotation phase <f>. The magnetic 
latitude of this vector is 6e- The vectors U and V are in the planes 
R-M and R-E respectively; for their application, see section [All 



i ^ 5 

4> — y uj 

Therefore the derivation of equation IA2I in section IA1I 
also demonstrates the correctness of equation [2] 

sin 0e — cos S sin Ob , A ^ 

coscj - — A4 

sin o cos 6 b 

In step |(vi)] we then wish to find the emission power P 
per solid angle Q in this direction. We may assume that the 
emission is evenly distributed around the cone (i.e. dP/duo 
is constant). If we differentiate equation IA4I with respect to 
0e, we get: 

= CQS ^ ( A5 ) 

cIQe cos Ob sin 5 sin uo 

This allows us to obtain the distribution of emission 
power with magnetic latitude: 

dP _ dP duj 
dOE duo dOE 

COS ^ E (as is constant) (A6) 



cos b sin S sin uj duo 
The emission power per solid angle is related to this: 



dP 1 1 dP 



dQ, 2n cos 0e dOE 
1 



cos Ob sin 5 sin 



(A7) 



This is the result used in step (vi) of section [3~4l equa- 
tion [3] 
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